Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 18 May 2010 



(MN style file v2.2) 



On the Size Distribution of Supernova Remnants in the 
Magellanic Clouds 

Carles Badenes^'^*, Dan Maoz^, Bruce T. Draine^ 

^School of Physics and Astronomy, Tel- Aviv University, Tel- Aviv 69978, Israel 

^ Benoziyo Center for Astrophysics, Weizmann Institure of Science, Rehovot 76100, Israel 

'"^Department of Astrophysical Sciences, Princeton University, Ivy Lane, Princeton, NJ 08540 



18 May 2010 



ABSTRACT 

The physical sizes of supernova remnants (SNRs) in a number of nearby galaxies fol- 
low an approximately linear cumulative distribution, contrary to what is expected for 
decelerating shock fronts. This phenomenon has been variously attributed to observa- 
tional selection effects, or to a majority of SNRs being in "free expansion" , with shocks 
propagating at a constant velocity into a tenuous ambient medium. We compile multi- 
wavelength observations of the 77 known SNRs in the Magellanic Clouds, and argue 
that they provide a fairly complete record of the SNe that have exploded over the last 
~ 20 kyr, with most of them now in the adiabatic, Sedov phase of their expansions. 
The roughly linear cumulative distribution of sizes (roughly uniform in a differential 
distribution) can be understood to result from the combination of the deceleration 
during this phase, a transition to a radiation-loss-dominatcd phase at a radius that 
depends on the local gas density, and a probability distribution of densities in the 
interstellar medium varying approximately as p~^. This explanation is supported by 
the observed powerlaw distributions, with index — 1, of three independent tracers of 
density: neutral hydrogen column density. Ha surface brightness, and star-formation 
rate based on resolved stellar populations. In this picture, the observed cutoff at a 
radius of 30 pc in the SNR size distribution is due to a minimum in the mean ambient 
gas density in the regions where supernovae (SNe) explode. We show that M33 has a 
SNR size distribution very similar to that of the Magellanic Clouds, suggesting these 
features, and their explanation, may be universal. In a companion paper (Maoz & 
Badenes 2010), we use our sample of SNRs as an effective "SN survey" to calculate 
the SN rate and delay time distribution in the Magellanic Clouds. The hypothesis that 
most SNRs arc in free expansion, rather than in the Sedov phase of their evolution, 
would result in SN rates that are in strong conflict with independent measurements, 
and with basic stellar evolution theory. 

Key words: supernovae: general - supernovae: remnants - galaxies: individual: LMC, 
SMC 



1 INTRODUCTION 

The ecology of galaxies is dominated by supernova (SN) 
explosions, which inject energy and enriched material into 
the interstellar medium and trigger the formation of the 
next generations of stars. Many fundamental aspects of SNe 
are still poorly understood, both for the core-collapse (CC) 
SN explosions that are thought to end the lives of massive 
(> 8Mq) stars and for the type-la SNe (SNe la) that are 
believed to be the thermonuclear combustions of CO white 
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dwarfs (WDs) that approach the Chandrasekhar mass. Nev- 
ertheless, the two flavors of SNe deposit a similar amount 
(~ 10^^ erg) of kinetic energy into the surrounding medium, 
leaving behind supernova remnants (SNRs) that remain vis- 
ible for thousands of years. In recent times, the study of 
young SNRs at X-ray wavelengths has emerged as a new 
way to e xplore the phys ics of CC and Type la SN explo- 
sions (see lBadenesll201(J . and references therein), but most 
known SNRs are too old to provide much information about 
th e specific events that originated them (se e discussions 
in iRakowski et all l2006l : iBadenes et"al] |2009| ). In spite of 
this, much can be learned by studying the entire population 
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of SNRs, young and old, in nearby galaxies. Because the 
timescales for SNR evolution are short compared to most of 
the processes that affect the structure of galaxies, SNR cata- 
logues provide a clean record of the environments where SNe 
explode, which can be used to put constraints on th e prop- 
erties of their progenitors (e.g. iBadenes et al.ll2009l ). More- 
over, by considering the properties of the entire population 
of SNRs in a galaxy together with the bulk properties of the 
gas they are expanding into, we can gain insights into the 
evolutionary phases of SNRs (|Woltieij|l972l ). the structure of 
galax ies on scales comparable to the average SNR size 
[2003), an d the cycles of matter and energy in the interstellar 
medium (|Ferrierj|20"0ll ). 

In this paper and in a companion publication (Maoz & 
Badenes 2010, henceforth Paper II), we use the SNR popula- 
tion in the Magellanic Clouds to explore some of these issues. 
The Magellanic Clouds (MCs) have the advantage of being 
close enough to study key aspects of their global structure 
in great detail, and they also harbor a large and extensively 
observed population of SNRs. Thus, they are the optimal 
setting to study the interplay between local density, star 
formation, SN explosions, and SNR evolution on a galactic 
scale. Our ultimate goal, and the focus of Paper II, is to 
derive the SN rate and delay time distribution (i.e., the SN 
rate as a function of time following a brief burst of star for- 
mation) in the Magellanic Clouds. However, this cannot be 
done without understanding first the relationship between 
the lifetime of SNRs and the properties of their local envi- 
ronments. This is the subject of the present work. 

The evolution of SN Rs has been the subject of many 



3e en tne suo.iect or many 

, ^ I1972I; fChevalieil [1982; 

Cioffi et al.l 1 19881 : iBlondin et al.l 1 19981 : iTruelove fc McKed 



theoretical studi e s (e.g.. IWoltje 



_ ■ _ ■ jTruelove , 

19991 ). Because accurate ages are only known for a handful 



of young, often historical objects, any observational tests of 
these theoretical models must rely on SNR size as a proxy 
for age. Given that SNRs of equal ages will have differ- 
ent sizes if they expand in different media, this necessar- 
ily brings the role of local density into the picture. Pre- 
vious works on the distribution of SNR sizes initially fo- 
cused on the Milky Way and the MCs (e .g. Mathewson et akl 
1984 iGreenI Il984|: iHughes et al.l 1 1984 iBerkhuiisenI 1 19871 : 



Chu fc Kennicut? 19881 ). but more recent efforts have also 



explored other galaxi es in the Local Group, inc luding M31 



Magnier et al.H 19971 ). M33 (|Long et all I2OIOI ). and M83 
Dopita et al.ll2010h . With few exceptions, these studies gave 



little consideration to the bulk properties of the gas in the 
galaxies hosting the SNRs. In many cases, their samples were 
also affected by issues of completeness and biases from work- 
ing at a single wavelength. Not surprisingly, these efforts 
have failed to produce a unified, physically motivated pic- 
ture of the evolution of SNRs in the interstellar medium. 
Here, we propose a first approximation to the problem in 
the context of the Magellanic Clouds. 

This paper is organized as follows. In § [51 we present 
a compilation of multi-wavelength observations for the 77 
known SNRs in the Magellanic Clouds, and we argue that 
it provides a fairly complete record of all the SNe that have 
exploded over the last ~ 20 kyr. In §[1 we examine the size 
distribution of SNRs in both galaxies, and we find that the 
cumulative distribution is close to linear (i.e., the differen- 
tial is close to uniform), within the uncertainties associated 
with the relatively small number of objects, up to a marked 



cutoff at a physical radius of ~ 30 pc. We also show that 
the SNR size distribution in M33 has very similar proper- 
ties, suggesting that these features might be widespread. In 
§|4l we propose a physical model to explain this distribution, 
based on the assumption that most objects are in the Sedov 
(adiabatic) stage of their evolution, and that they rapidly 
fade away once they transition to the radiative stage, at an 
age that depends on the local density. Under these condi- 
tions, the uniform distribution of SNR sizes requires that 
the gas density in the Clouds have a probability distribu- 
tion described by a power law with an index of —1 (i.e., 
5P/5p ~ P~^). In § [S] we test this requirement by exam- 
ining the distribution of three independent density tracers 
in the Clouds: neutral hydrogen column density. Ho? surface 
brightness, and star-formation rate based on resolved stellar 
populations. We find that these tracers are indeed well de- 
scribed by powerlaws with a —1 index. This lends credence 
to our model, and provides us with the crucial means to 
estimate the visibility time of SNRs in different locations, 
which we review briefly in § |6l and more extensively in Pa- 
per II, where we use it to derive the SN rate and delay time 
distribution in the Clouds. We conclude by summarizing our 
main results and outlining avenues for future work in § [T] 



2 THE SUPERNOVA REMNANTS IN THE 
MAGELLANIC CLOUDS: OBSERVATIONS 
AND SAMPLE COMPLETENESS 

The population of SNRs in the MCs has been the object of 
extensive study for many decades. Several catalogues have 
been compiled at difi'erent wavelengths, from the radio to 
the optical and X-ray. Because SNRs in the Milky Way and 
the MCs are usually discovered in the radio, it is the ra- 
dio catalogues that often have the largest number of entries, 
but there is some confusion in the li terature regarding the 
total number of SNRs in the Clouds. iFilipovic et al.l (| 19981 ) 
listed all the discrete radio sources in the Parkes survey of 
the MCs, and found 32 SNRs and 12 SNR cand idates in the 
LMC, and 12 SNRs in the SMC. According toU^neeTal] 
([2003), a revision of the Parkes survey complemented with 
Australian Telecope Compact Array (ATCA) data yields 52 
confirmed SNRs and 20 candidates in the LMC, but these 
sources are not listed in their paper. Instead, the authors 
present optical spectroscopy of 25 of the 52 confirmed LMC 
SNRs. Data from ATCA wa s also collected for the SMC, 
where IFilipovic et al.l (|2005l ) list 16 confirmed SNRs. An 
optical catalogue of MC SNRs is being assembled as part 
of the Magellanic Clouds Emission Line Survey (MCELS, 
ISmith et al1l2000l ): the most recent on-line version of this 
catalogutEl fist s 31 SNRs in the LM C and 11 in the SMC. 
In the X-rays, IWiUiams etHI |l99i) published an atlas of 
ROSAT sources, which contained 31 LMC SNRs. The most 
recent on-line version of this catalogu^l, complemented with 
optical, radio and in frared data, lists 38 c o nfirm ed SNRs 
in the LMC. Finally, Ivan der Hevden et al.l (|2004l ) present 
XMM-Newton observations of 13 SMC SNRs. 



^ Last modified January 9 2006, see 
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Figure 1. Maps of the LMC (top, from lKim et al.ll2003h and the SMC (bottom, from IStanimirovic et"allll999h in the HI 21 cm Une. 
The positions of the SNRs from Table 1 are indicated by crosses: 54 objects in the LMC, 23 objects in the SMC. The contours are at 
column densities of 5, 10, 30, 50, and 100 X 10^0 cm-^. 
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Taken separately, these catalogues are not adequate 
for our purposes, because they disagree on the classifica- 
tion of particular sources, and assign different names, sizes 
and positions to well-known objects. Catalogues at differ- 
ent wavelengths can also be sensitive to SNRs at different 
stages of their evolution, so a multi-wavelength compilation 
is the best way to obtain a global picture of the SNR pop- 
ulation. At the writing of this paper, there is still no uni- 
fied, multi-wavelength catalogue of all the SNRs in the MCs 
in the refereed literature, although the Magellanic Clouds 
Sup ernova Remnants (MCSNR) c ollaboration is assembling 
one (|Murphv Williams et al.ll2010l ) 0. The current version of 
the MCSNR catalogue contains 48 confirmed objects in the 
LMC and f9 in the SMC. We have merged this on-line cata- 
logue with the ones listed in the previous paragraph, taking 
only the sources that are classified as confirmed SNRs in 
at least one catalogue, and eliminating all duplicates and 
sources listed only as candidate SNRs. This yields 54 con- 
firmed SNRs in the LMC and 23 in the SMC, which we list 
in Table 1. In Figure [1] we plot the positions of these 77 
SNRs on t he HI maps of the MCs p ublished bv iKim et al.l 
l|2003h and lStanimirovic et al.l (| 19991 '). 

Because we are interested mostly in reliable sizes (to 
derive the histograms in § [3} and positions (to calculate 
the delay time distribution in Paper II), we take the cen- 
tres and diameters of the SNRs from high resolution X-ray 
observations {Chandra or XMM-Newton), whenever possi- 
ble. If no high-resolution X-ray data are available, we revert 
to the sizes and centers listed in the original catalogues, 
which were determined on a case- by-case basis, usually se- 
lecting the wavelength that provided the cleanest measure- 
ment (see Table 1 and references for details). We have found 
that high-resolution X-ray data provide the most reliable 
measurements for several X-ray bright and radio faint ob- 
jects, which are often small and have boundaries that are 
hard to resolve clearly using single-dish radio data. For the 
bulk population, sizes determined at different wavelengths 
might disagree in individual cases, but this does not have a 
large impact on the final distribution of sizes, as long as the 
data sets are of rea sonable quality (e.g., see Figure 7 and 
related discussion in iFilipovic et aLlbOOSl ). 

With only one exception (SNR J0051. 9-7310 in the 
SMC, see Table 1), all the SNRs in the MCs are de- 
tected in the radio, and their flux densities can be found 
either in the MCSNR catalogu e or in the general m ulti- 
band radio source catalogues of IFilipovic et all (| 19951 ) and 
IFilipovic et"all (|2002l '). We have listed these flux densities in 
Table 1, at 1.4 GHz when available, and at higher frequen- 
cies when there was no 1.4 GHz detection. In Figure (21 we 
plot the radio flux densities as a function of SNR dia meter, 
adop ting the flducial distances of 50 kpc to the LM C (jAlvesI 
|2004| ) and 60 kpc to the SMC (|Hilditch et al.ll2005^ to con- 
vert angular sizes to linear dimensions (we will use these 
distances for the remainder of the paper). The flux densities 
of the SMC SNRs have been multiplied by a factor 1.44 for 
comparison to the LMC SNRs. Remarkably, there appears 
to be a "floor" in SNR flux density at ~ 50 mjy, with only 
a few objects being fainter than this limit. 



^ For the most recent on-line version (2009), see 
|http : //hoth. ccssc ■ org/mcsnrTI 
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Figure 2. Integrated radio flux densities of SNRs in the LMC 
(empty symbols) and the SMC (filled symbols), plotted versus 
SNR diameter. The flux densities of the SMC SNRs have been 
multiplied by 1.44 for comparison. The circles indicate flux densi- 
ties at 1.4 GHz, the triangles are at higher frequencies (see Table 
1 for details). The dashed and dotted lines show conservative So" 
limits for the selection of candidate SNRs in both galaxies, ac- 
counting for the fiux limit per beam and the number of beams per 
objec t of given diameter (see iHughes et al.|[2007l : IFilipovic et al.l 
l2002h . The SMC limit has also been scaled by 1.44. There is a 
real paucity of SNRs below a flux density of ~ 50 mJy at 50 kpc, 
corresponding to a luminosity of ~ 1.5 X 10^^ W Hz~^. 



The degree of completeness of our sample is hard to 
estimate, because we include objects that were found using 
a very hetrogeneous mix of selection criteria, but a num- 
ber of factors indicate that we should not be missing a 
large number of SNRs. First, the floor at 50 mJy is much 
brighter than the nominal sensitivity of the ATCA/Parkes 
radio surveys, which at 1.4 GHz is around 0.3 mJy per beam 
in the LMC jHuehes et al.l|2007l') . and 1.8 mJy per beam 

20021 ). with beam sizes of 40" 
respectively. Even requiring at least a 5a" detec- 



m the SMC (|Filipovic et al 

and 98' 



tion for the selection of SNR candidates, and allowing for 
the fact that large SNRs should be harder to identify at 
the same level of significance, the vast majority of our ob- 
jects are safely above any conservative detection limits in 
the radio (see Figure [2]). Second, the syst ematic searches 
for SN Rs o n the ATCA / Parke s surveys bv IFilipovic et al] 
(j2005l ) and IPavne et al.l ()2008l ) have not uncovered a sub- 
stantial population of previously unknown SNRs between 
the radio floor and the detection limits. We note that the 
tot al number o f confi rmed SNRs reported (but not listed) 
bv lPavne et al.l (|2008l ) in the LMC is very close to our own 
count (52 vs. 54). Third, the existence of a substantial pop- 
ulation of radio detected SNRs that are masquerading as 
other kinds of sources is unlikely in the presence high quality 
multi-wavelength data. Small, radio-faint SNRs that might 
be mistaken for background sources are often young ob- 
jects with bright X-ray emission like SNRs J0509.5— 6731 
(B0509-67.5) or J0505. 7-6753 (DEM L71). Older, more 
extended objects that might be confused with surround- 
ing or nearby HII regions are usually identified at optical 
wavelengths by diagnostic quantities like the [S II] /Ha ra- 
tio, which can be used to distinguis h shocked nebulae from 
photoionized gas (|Fesen et al.lll985l ). 

A related concern is that some SNRs might be "lost" 
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Table 1. SNRs in the Magellanic Clouds 



SNR« 




Alternative 




Position (J2000) 




Diameter 


Wave- 


Refe- 


Radio flux 


Catalogue 






Name^ 




R.A. 




Decl. 




(arcsec) ^ 


length^ 


rence^ 


density (Jyi£ 


entries^ 


LMC SNRs 


J0448.4- 


-6660 




04h 


48m 


22s 


-66d 


59m 


52s 


220 


R 


M 


0.04 


M 


J0449.3- 


-6920 




04h 


49m 


20s 


-69d 


20m 


20s 


133 


R 


M 


0.11 


M 


J0450.2- 


-6922 


B0450-6927 


04h 


50m 


15s 


-69d 


22m 


12s 


210 


R 


P 


0.27^^ 


P 


J0450.4- 


-7050 


B0450-709 


04h 


50m 


27s 


-70d 


50m 


15s 


357 


X 


W 


o.se" 


WSFM 


J0453.2- 


-6655 


N4 


04h 


53m 


14s 


-66d 


55m 


13s 


252 


X 


W 


0.11 


WPFM 


J0453.6- 


-6829 


B0453-685 


04h 


53m 


38s 


-68d 


29m 


27s 


120 


X 


G03 


0.19 


WFSM 


J0453.9- 


-7000 


B0454-7005 


04h 


53m 


52s 


-70d 


00m 


13s 


420 


R 


P 


0.262. 


P 


J0454.6- 


-6713 


N9 


04h 


54m 


33s 


-67d 


13m 


13s 


177 


X 


S06 


o.o" 


WFM 


J0454.8- 


-6626 


NllL 


04h 


54m 


49s 


-66d 


25m 


32s 


87 


X 


W 


0.11 


WPFSM 


J0455.6- 


-6839 


N86 


04h 


55m 


37s 


-68d 


38m 


47s 


348 


X 


W 


0.26 


WPFSM 


J0459.9- 


-7008 


N186D 


04h 


59m 


55s 


-70d 


07m 


52s 


150 


O 


W 


0.07 


WPFSM 


J0505.7- 


-6753 


DEM L71 


05h 


05m 


42s 


-67d 


52m 


39s 


72 


X 


B07 


0.01 


WSM 


J0505.9- 


-6802 


N23 


05h 


05m 


55s 


-68d 


01m 


47s 


111 


X 


H06 


0.35 


WPFSM 


J0506.1- 


-6541 




05h 


06m 


05s 


-65d 


41m 


08s 


408 


R 


M 


0.11 


M 


J0506.8- 


-7026 


B0507-7029 


05h 


06m 


50s 


-70d 


25m 


53s 


330 


R 


P 


0.36£ 


P 


J0509.0- 


-6844 


N103B_^ 


05h 


08m 


59s 


-68d 


43m 


35s 


28 


X 


B07 


0.51 


WPFSM 


J0509.5- 


-6731 


B0509~67.5 


05h 


09m 


31s 


-67d 


31m 


17s 


29 


X 


B07 


0.08 


WSM 


J0513.2- 


-6912 


DEM L109 


05h 


13m 


14s 


-69d 


12m 


20s 


215 


R 


Boj07 


0.20 


WFSM 


J0518.7- 


-6939 


N120 


05h 


18m 


41s 


-69d 


39m 


12s 


134 


X 


R08 


0.35 


WFM 


J0519.6- 


-6902 


B0519-690 


05h 


19m 


35s 


-69d 


02m 


09s 


31 


X 


B07 


0.10 


WFSM 


J0519.7- 


-6926 


B0520-694 


05h 


19m 


44s 


-69d 


26m 


08s 


174 


O 


W 


0.12 


WPFSM 


Tn'i2i fi- 


-6543 




05h 


21m 


39s 


— 65d 


43m 


07s 


I 


R 


M 


0.03 


M 


Tn'i23 1 - 


-6753 


N44 


05h 


23m 


07s 


— 67d 




12s 


228 


R 


K98 


0.14 


WPFM 

V V i- i- 1V± 




-6624 


DEM T,1 75a 


05h 


24m 


20s 


— 66d 


24m 


23s 


234 


o 


W 


0.11 


WPFSM 

V V J- J- 1V± 




-6938 


N132D 


05h 


25m 


04s 


— 69d 


38m 


24s 


114 


o 


W 


3.71 


WPFSM 

V V J- J- kj 1V± 


70525 4- 


-6559 


N49B 


05h 


25m 


25s 


— 65d 


59m 


19s 


168 


X 


P03b 


0.32 


WPFSM 




-6605 


N49 


05h 


26m 


00s 


— 66d 


04m 


57s 


84 


X 


P03 


1.19 


WPFSM 

V V J- J- tJ i v± 


7(1527 fi- 


-6912 




05h 


27m 


39s 


— 69d 


12m 


04s 


147 


o 


W 


0.05 


WPFSM 

V V J- J- kJ 1V± 


J0527.9- 


-6550 


DEM L204 


05h 


27m 


54s 


-65d 


49m 


38s 


303 


o 


w 


0.14 


WPFSM 


J0527.9- 


-6714 


B0528-6716 


05h 


27m 


56s 


-67d 


13m 


40s 


196 


R 


M 


0.10 


FSM 


J0528.1- 


-7038 


B0528-7038 


05h 


28m 


03s 


-70d 


37m 


40s 


60 


R 


P 


0.28.^ 


P 


J0529.1- 


-6833 


DEM L203 


05h 


29m 


05s 


-68d 


32m 


30s 


667 


R 


M 


0.24 


M 


J0529.9- 


-6701 


DEM L214 


05h 


29m 


51s 


-67d 


01m 


05s 


100 


R 


P 


0.06 


PM 


J0530.7- 


-7008 


DEM L218 


05h 


30m 


40s 


-70d 


07m 


30s 


213 


R 


M 


0.06 


M 


J0531.9- 


-7100 


N206 


05h 


31m 


56s 


-71d 


00m 


19s 


192 


O 


W 


0.33 


WFSM 


J0532.5- 


-6732 


B0532-675 


05h 


32m 


30s 


-67d 


31m 


33s 


252 


X 


W 


0.16 


WSM 


J0534.0- 


-6955 


B0534-699 


05h 


34m 


02s 


-69d 


55m 


03s 


114 


X 


H03 


0.08 


WPFSM 


J0534.3- 


-7033 


DEM L238 


05h 


34m 


18s 


-70d 


33m 


26s 


180 


X 


B06 


0.06 


WPFSM 


J0535.5- 


-6916 


SNR1987A 


05h 


35m 


28s 


-69d 


16m 


lis 


2 


R 


N08 


0.05 


WM 


J0535.7- 


-6602 


N63A 


05h 


35m 


44s 


-66d 


02m 


14s 


66 


X 


W03 


1.43 


WPFSM 


J0535.8- 


-6918 


Honeyeomb 


05h 


35m 


46s 


-69d 


18m 


02s 


102 


X 


W 


0.30 


WFSM 


J0536.1- 


-6735 


DEM L241 


05h 


36m 


03s 


-67d 


34m 


36s 


135 


X 


Ba06 


0.29 


WPFSM 


J0536.1- 


-7039 


DEM L249 


05h 


36m 


07s 


-70d 


38m 


37s 


180 


X 


B06 


0.05 


WPFSM 


J0536.2- 


-6912 


B0536-6914 


05h 


36m 


09s 


-69d 


11m 


53s 


480 


R 


P 


4.26" 


PF 



" SNR names are constructed from the location of the centre in J2000 coordinates as given in the catalogues. 

" Alternative names are given for clarification purposes; see IWilliams et al. I lll999h and iFilipovic et al. for a discussion 

of past conventions in SNR nomenclature. 

Some catalogues give two axial diameters for asymmetric SNRs. In those cases, we list the mean diameter. 

For the center location and diameter. R: radio; O: optical; X: X-ray 
^ For the center location and diameter. See table notes for a list of all the catalogues and other references. 
f All flux densities are at 1.4 GHz from catalogue M, unless indicated. 
3 Catalogues that list the object as a confirmed SNR; see table notes. 
^ At 4.75 GHz, from F95 
* At 2.45 GHz, from F95 
3 From F95 

^ This widely used name is the result of an incorrect identification, see entry in M 
^ No size is given for this SNR in the MCSNR catalogue. 

At 2.45 GHz, from F95 
" From F95 
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Table 1 - continued 



LMC SNRs 



J0537.4-- 


^6628 


DEM L256 


05h 


37m 27s 


-66d 


27m 


50s 


204 


R 


M 


0.06 


M 


J0537.6- 


-6920 


B0538-6922 


05h 


37m 37s 


-69d 


20m 


23s 


169 


O 


S 


1.01^ 


FS 


J0537.8- 


-6910 


N157B 


05h 


37m 46s 


-69d 


10m 


28s 


102 


X 


C06 


2.64" 


WFSM 


J0540.0- 


-6944 


N159 


05h 


39m 59s 


-69d 


44m 


02s 


78 


X 


WOO 


1.90 


WM 


J0540.2- 


-6920 


B0540-693 


05h 


40m lis 


-69d 


19m 


55s 


60 


X 


HOI 


0.88 


WFSM 


J0543.1- 


-6858 


DEM L299 


05h 


43m 08s 


-68d 


58m 


18s 


318 


X 


W 


0.21 


WSM 


J0547.0- 


-6943 


DEM L316B 


05h 


46m 59s 


-69d 


42m 


50s 


84 


X 


W05 


0.52 


WFSM 


J0547.4- 


-6941 


DEM L316A 


05h 


47m 22s 


-69d 


41m 


26s 


56 


X 


W05 


0.33 


WFM 


J0547.8- 


-7025 


B0548-704 


05h 


47m 49s 


-70d 


24m 


54s 


102 


X 


H03 


0.05 


WFSM 


J0550.5- 


-6823 




05h 


50m 30s 


-68d 


22m 


40s 


312 


R 


M 


0.64 


M 


SMC SNRs 


J0040.9- 


-7337 


DEM S5 


OOh 


40m 


55s 


-73d 


36m 


55s 


121 


R 


M 


0.19 


M 


J0046.6- 


-7309 


DEM S32 


OOh 


46m 


39s 


-73d 


08m 


39s 


136 


X 


H 


0.07 


HMF05 


J0047.2- 


-7308 


IKT2 


OOh 


47m 


12s 


-73d 


08m 


26s 


66 


X 


H 


0.46 


HMF05 


J0047.5- 


-7306 


B0045-733 


OOh 


47m 


29s 


-73d 


06m 


01s 


180 


R 


F05 


0.14^ 


F05 


J0047.7- 


-7310 


HFPK419 


OOh 


47m 


41s 


-73d 


09m 


30s 


90 


X 


H 


O.l" 


HM 


J0047.8- 


-7317 


NS21 


OOh 


47m 


48s 


-73d 


17m 


27s 


76 


R 


F05 


0.03 


F05 


J0048.1-^ 


-7309 


NS19 


OOh 


48m 


06s 


-73d 


08m 


43s 


79 


R 


F05 


0.08 


F05 


J0048.4- 


-7319 


IKT4 


OOh 


48m 


25s 


-73d 


19m 


24s 


84 


X 


H 


0.15 


HMF05 


J0049.1- 


-7314 


IKT5 


OOh 


49m 


07s 


-73d 


14m 


05s 


116 


X 


H 


0.28 


HMF05 


J0051.1- 


-7321 


IKT6 


OOh 


51m 


07s 


-73d 


21m 


26s 


144 


X 


H05 


0.10 


HMF05 


J0051.9- 


-7310 


IKT7 


OOh 


51m 


54s 


-73d 


10m 


24s 


97 


X 


M 


0.00^ 


M 


J0052.6- 


-7238 


B0050-728 


OOh 


52m 


33s 


-72d 


37m 


35s 


323 


R 


F05 


0.22" 


MF05 


J0058.3- 


-7218 


IKT16 


OOh 


58m 


16s 


-72d 


18m 


05s 


200 


X 


H 


0.07 


HMF05 


J0059.4- 


-7210 


IKT18 


OOh 


59m 


25s 


-72d 


10m 


10s 


158 


X 


H 


0.37 


HMF05 


J0100.3- 


-7134 


DEM S108 


Olh 


00m 


21s 


-71h 


33m 


40s 


149 


R 


F05 


0.21 


MF05 


J0103.2- 


-7209 


IKT21 


Olh 


03m 


13s 


-72d 


08m 


59s 


62 


X 


H 


0.10 


HMF05 


J0103.5- 


-7247 


HFPK334 


Olh 


03m 


30s 


-72d 


47m 


20s 


86 


R 


M 


0.05 


M 


J0104.0- 


-7202 


B0102-7219 


Olh 


04m 


02s 


-72d 


01m 


48s 


44 


X 


GOO 


0.28 


HMF05 


J0105.1- 


-7223 


IKT23 


Olh 


05m 


04s 


-72d 


22m 


56s 


170 


X 


P03c 


0.09 


HMF05 


J0105.4- 


-7209^ 


DEM S128 


Olh 


05m 


23s 


-72d 


09m 


26s 


124 


X 


H 


0.05 


HMF05 


J0105.6- 


-7204"" 


DEM S130 


Olh 


05m 


39s 


-72d 


03m 


41s 


79 


R 


F05 


0.05 


F05 


J0106.2- 


-7205 


IKT25 


Olh 


06m 


14s 


-72d 


05m 


18s 


110 


X 


H 


0.01 


HMF05 


J0114.0- 


-7317 


N83C 


Olh 


14m 


00s 


-73d 


17m 


08s 


17 


R 


M 


0.02'= 


M 



C atalogues: 

F: lFilipovic et ahl lll998f). Catalogue of discrete radio sources in the MCs. Data from Parkes telescope. 

F95: ^^lipovlc et al. I lll995l '). Catalogue of discrete radio sources in the LMC, with fluxes at 1.40, 2.45, 4.75, 4.85 and 8.55 GHz. 
Data from Parkes te l escop e. 

F02: lFilipovic et al.l ||2002| '). Catalogue of discrete radio sources in the SMC, with fluxes at 1.42, 2.37, 4.80, and 8.64 GHz. Data 
from ATCA and Parkes tel escope. 

F0 5:lFilipovic et al.l bOO^). Ca talogue of radio SNRs in the SMC. Data from ATCA. 

H: lvan der Hevden et al.l ||2004|) . Catalogue of X-ray bright SNRs in the SMC. Data from XMM-Newton. 
M: MCSNR on-line d atabase, (R. Williams et al.) [http : //both. ccssc . org/mcsnr/ 

P: iPavne et al.l 1 I2OO8I) . Optical spectroscopy of radio SNRs in the LMC. Radio data from ATCA, optical data from Siding 
Spring Observatory and South African Astronomical Observatory. 

S: MCELS on-line ca talogu e (C. Smith et al.), http : //www, cti o .noao . edu/mcels/snrs/snrcat .html | Optical data from CTIO. 
W: IWilliams e t~ar ('1999"). X-ray atlas from ROSAT, also including radio and optical observations. Available on-line at 
|littp: //www. astro . Illinois . edu/projects/atlas/index .html 

Other Referen ces: 

B06: ,Borkowski et al.l 1 I2OO6I): Ba06: iBamba et al.l ll2006h: B07: iBadenes et al.l ||2007^: Bo.i0 7: iBoiicic et al.l ll2007t); C06 



IChen et al.l bOOell: GOO 



Gaetz 



IWilliams et al.l 1I2OOOI '): W03: IWarren et al.l | |2003|) : W05: IWilliams fc Chul ||2005|) " 



~ ^''^ llyOOl^'GOS: ICaensler et all ||2003|'): HOI: iHwane et al.l bOOll); H0 3: iHendrick et al, 

20031): H05: IHendrick et al.l ||2005^ : HOe TlHughes et a l.l (|2006|): K" 98:lKim et al.l lll998h: N08:lNg et'alTll2008t): P03:|Park etaL. 
2003al'): P03b:|Park c t al. (20 03d'): P03c:|Park et al.l (2003b'): R08:lReves-Iturbide et al.l ll2008i' ): S06: ISeward et al] 1I2OO6I) : WOO: 



At 4.75 GHz, from F95 
From F02 

"= This SNR has no radio or optical counterpart 

This source might be two smaller SNRs, sec H, F05 
^ At 50 GHz 
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inside superbubbles formed by several nested SN explosions 
(IMac Low fc McCravHTosi '). In practice, SNRs evolving in- 
side low-density cavities do not disappear - they just be- 
come fainter and evolve more slowly, at least until the blast 
wave reaches the cavity edges. The largest, most prominent 
superbubble complex in the LMC, 30 Doradus, harbors sev- 
eral well-studied SNR s, including J0537.8-6910 (N157B) 
iTownslev et al.ll2006h . and many other objects listed in Ta- 
ble 1 are associated with superbubbles. As long as the sur- 
veys are sensitive enough to find faint SNRs, this should not 
be a major issue. 

In the absence of a systematic re-analysis of all the 
available data with a set of homogeneous criteria to identify 
SNRs in both Clouds, we cannot guarantee that our sam- 
ple is 100% complete. Such an analysis is outside the scope 
of the present work, but all the evidence seems to indicate 
that our compilation should not be missing a large number 
of objects, and it should include few, if any, spurious ones. 
We conclude that our list of MC SNRs is at least fairly com- 
plete (which is enough for our goals), and that the paucity 
of SNRs below ~ 50 mjy is probably real, and not due to 
any observational selection effects. This suggests that any 
MC SNRs with flux densities between the radio floor and 
the detection limits must fade relatively quickly. As we shall 
discuss in § [51 the results of Paper II can be used to argue 
independently that the SNR sample we have compiled has 
a high degree of completeness. 



3 THE DISTRIBUTION OF REMNANT SIZES 

Figure [3] shows the cumulative and differential histograms of 
remnant sizes in the LMC and in the SMC. In both Clouds, 
the cumulative distributions are roughly linear in remnant 
size, i.e., with an equal number of remnants per equal size 
bin in the differential distributions, up to a cutoff at a phys- 
ical radius rcut ~ 30 pc. 

Because the data points in a cumulative histogram are 
not independent of one another, these curves cannot be fit- 
ted in the usual way (using the statistic). A robust, quan- 
titative estimate for the slope of the distributions can be 
obtained by pe rforming instead a maximum-likelihood fit, 
as follows (see iMaoz fc Rixl 1 19931 . for a similar treatment 
applied to a different problem). Suppose that a particular 
model predicts a size distribution dN/dr = n{r), which in- 
tegrates to /J^""' n{r)dr = N, where A'' is the total number 
of remnants up to rcut. If we bin our data into many small 
bins between r = and rcut, each of width Sr, most bins 
will have zero SNRs, and some will have one SNR. Given 
the model, the Poisson probability of finding j remnants in 
the ith bin, for which the model predicts n{ri) remnants, is 



PUHn)) 



-7i(r j ) 



n{riY/jl 



(1) 



We will define the likelihood of a given model as the prod- 
uct of these probabilities. The logarithm of the likelihood, 
considering that j always equals either or 1, simplifies to 



lnL = [ ^ \nn{ri)]-N, 

i(j>0) 



(2) 



where the summation is only over the specific data values of 
the SNR radii. A power law of index a, having the proper 



normalization, will have the form 

, , Nja + l) , 



(3) 



Inserting in Eq.[2l differentiating InL with respect to a, and 
equating to zero to find the maximum gives the maximum 
likelihood solution. 



(q + l)ml = 



N 



A In rcut - Ei>oln 



with an uncertainty on a of 



Aa = 



d^(lnZ,) 
da2 



-1/2 



a + 1 



N 



(4) 



(5) 



This procedure yields a maximum likelihood index of 
a = 0.14 ± 0.18 for the LMC, and a = 0.32 ± 0.28 for 
the SMC. Thus, the LMC appears to indeed have a SNR 
size distribution that is close to uniform. In the SMC, the 
best fit is intermediate between a fiat distribution and one 
that rises linearly with radius, but given the smaller number 
of SNRs, it is consistent with both slopes. From Figure [3l 
it appears that the steeper slope in the SMC is driven by 
the small-radius side of the distribution. Indeed, if we fit 
separately the first 8 points and the following 14 points, the 
maximum likelihood solution is a = 1.7 ± 1.0 at small radii, 
and a — 0.17 ± 0.23 thereafter. This result confirms the 
visual impression, but formally it is still consistent with a 
slope close to zero at all radii at the 1.7a level. We conclude 
that the SNR size distribution in both Clouds is consistent 
with being roughly uniform, although there are hints for a 
deviation at small radii in the SMC. It is possible that a 
deficit at small radii is also present in the LMC distribution 
(see Figure [3|, but the Poisson errors are too large to claim 
that the data require it. 

The linear cumulative distribution of SNR sizes 
in the MCs (and also the Mil ky Way) has bee n prev i- 
ously note d and discussed bv [ Mathewson et al.l ( 1984 ) , 
Mills et al.l h984). [Creenl (Il984l ). iHughes et aTl msj). 
Fusco-Femiano fc P roitc-Martincj ( 1984h . iBerkhuiisenI 
(Il987f). IChu fc Kennicutt (1988,), and most recently by 
iBandiera fc Petruk ( 20101 '). Several of these papers also 



pointed out cutoffs in the distribution. All of these papers 
considered smaller SNR samples, often based on much 
shallower radio data and, with few exceptions, did not 
include multi-wavelength observations. Some of these au- 
thors interpreted the observed size distribution as evidence 
that most MC SNRs are in their "free expansion" phase, 
during which the shock velocity is constant, and inferred 
that these SNRs expan d into an ex tremel y low-density 
medium. Alternatively, ICreenl (| 19841 ') and iHughes et al.l 
(|l984h warned that the observed distribution was the result 
of selection effects; most objects they discussed had been 
selected in X-rays, and the X-ray flux limits then led to 
the exclusion of larger and fainter remnants, and their 
faint radio counterparts. Our present compilation is bigger, 
it incorporates the most recent multi-wavelength data, it 
extends to larger sizes, and most importantly, it is sensitive 
to radio flux densities two orders of magnitude below the 
observed luminosity floor. With these data we now confirm 
the luminosity fioor, the uniform size distribution, and the 
cutoff at rcut ~ 30 pc in the Magellanic Cloud SNRs. 

These features of the SNR size distribution are also 



8 Badenes, Maoz, & Draine 



1 — ' — ' — ' — I — ' — ' — ' — I — ' — ' — "Tn — ' — ' — ' — I — ' — ' — ' — q -| 2 I — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — r 




Diameter [pc] Diameter [pc] 

Figure 3. Distribution of the SNR sizes listed in Table 1 for the LMC (red solid plots) and SMC (blue dashed plots), represented as 
cumulative (left) and differential (right) histograms. For illustrative purposes, linear size distributions normalized to give the correct 
number of SNRs at a diameter of 40 pc have been overplotted as dotted lines on the cumulative histograms. Both galaxies have SNR 
size distributions that are close to linear in the cumulative, or flat in the differential, up to a cutoff at r ~ 30 pc. 



present in other galaxies. Due to its proximity and face-on 
orientation, M33 has probably the best SNR sample outside 
of the MCs. The most recent catalo gue of M33 SNRs has 
been published bv lLong et al.l (|2010l ). and it includes data 
in the radio, opt ical, and X-rays (from the ChASeM33 sur- 
vey by Chandra, IPlucinskv" et al.ll2008l). The dis tribution of 
SNR sizes given in Table 3 of lLong et aD (j2010l ). which we 
plot in Figure [4] has the same broad characteristics we find 
in the Clouds, but backed by much better statistics (137 ob- 
jects). The cutoff at a physical radius of r ~ 30 pc is very 
clear, and the distribution appears remarkably uniform for 
radii between ~ 10 pc and the cutoff. In this size range, our 
method to estimate the maximum likelihood index yields 
a = 0.04 ± 0.11. Graphically, a linear function with the 
correct normalization reproduces the cumulative histogram 
quite well (see Figure |4| . A significant deficit of SNRs at 
radii below 10 pc is also obvious in the data. iLong et al.l 
(|2010l ) warn that a few small (young) SNRs might have es- 
caped detection (see their Section 9), but given the numbers 
involved, most of the deficit may be real. This indicates a 
deviation at small diameters from the uniform distribution 
along the lines of the one we noted in the MCs. 



4 PHYSICS BEHIND THE SIZE 
DISTRIBUTION 

The evolution of SN Rs is usually separated into three phases 
(e.g. IWoltieJll972l ). In the first, free-expansion, phase, the 
mass swept up by the expanding shock is small compared to 
the ejected mass, and the shock velocity is constant with 
time. In the Sedov-Tayloi0 phase, the swept-up mass is 
larger than the ejecta mass, and the expansion decelerates. 
However, the cooling time of the shocked gas is still longer 
than the age of the remnant, and the evolution is approx- 
imately adiabatic. Once the cooling time becomes compa- 
rable to the age, the SNR quickly loses energy radiatively, 
entering the radiation-loss-dominated snowplough phase, af- 
ter which it slows down, breaks up, and merges with the 
interstellar medium. 

Quantitative estimates show that typical SNRs should 
be in their Sedov- Taylor phases for ages between a few and 
a few tens of kyrs, and for sizes of order a few to a few 
tens of parsecs (e.g. ICioffi et al.lll988l : iBIondin et aLlll998l : 



iTruelove fc McKedll999l ). Thus, most of the MC SNRs in 



our sample should be in the Sedov- Taylor phase of their 
expansions, with radii growing as 



pl/5 -1/5,2/5 



(6) 



The available data seem to indicate that the distribu- 
tion of SNR sizes in M33 and the Magellanic Clouds is in- 
deed close to uniform between r ~ 10 pc and a sharp cutoff 
at Tcut ~ 30 pc. In the case of M33, this conclusion needs 
to be validated by a careful assessment of the completeness 
of the SNR s a mple, although from the discussion in § 7.3 
of iLong et all (|2010i ) it seems unlikely that a large number 
of M33 SNRs have escaped detection. Unfortunately, SNR 
samples of this quality are hard to obtain for other galax- 
ies, so it is difficult to verify how widespread these features 
of the SNR size distribution are in reality. We will return 
briefly to this issue in § [T] In the following section, we ar- 
gue that a uniform size distribution, including a deficit of 
SNRs at small radii and a sharp cutoff at a certain radius, 
arises naturally from the physics of SNR evolution and the 
distribution of densities in the interstellar medium. 



where Eo is the kinetic energy of the explosion, p is the 
ambient gas density, and t is the time. The shock velocity 
therefore decreases as 



V- t.^ p t , 

or equivalently expressed in terms of r rather than t, 

j-nl/2 -1/2 -3/2 



(7) 



(8) 



Reasonably assuming a constant SN rate in the LMC over 



^ The well-known solution for the propagation of a blast wave 
following a point-like explosion, commonly attributed to Leonid 
Sedov and Geoffrey Taylor, was also derived independently by 
John von Neumann. 
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Diameter [pc] Diameter [pc] 

Figure 4. Distribution of SNR sizes in M33 from I Long et al. represented as cumulative (left) and differential (right) histograms. 

A linear size distribution for diameters between 20 and 60 pc, normalized to give the correct number of SNRs at a diameter of 40 pc, 
has been overplotted as a dashed line on the cumulative histogram. 



the past few kyrs, dN/dt — const., the naively expected size 
distribution of SNRs in their Sedov phase is 



dN 
dr 



dN dt 
dt dr 



3/2 



(9) 



where we have used Eq. |8] to substitute for dt/dr. This is 
contrary to the observed nearly uniform [dN /dr ~ r'^) dis- 
tribution. As noted above, this has been interpreted before 
either as evidence for free expansion, or for selection effects, 
where fading of the older (and hence larger) remnants takes 
them below the detection limits, and culls the pileup in the 
distribution expected at large sizes for a decelerating pop- 
ulation. Both explanations are unlikely for our MC SNR 
sample. 

We propose, instead, that the uniform size dis- 
tribution arises as a result of the transition from 
the Sedov phase to the radiative phase (see also 
Fusco-Femiano fc Preite-Martinej|l984l : [Bandiera fc PetrukI 
The radius of this transition depends on ambient den- 
sity, with SNRs evolving in dense environments becoming 
radiative sooner than SNRs in low-density regions. When 
coupled with the distribution of densities in the Clouds, this 
leads to fewer and fewer sites at which large-radius Sedov- 
phase SNRs can exist. We show this as follows. 

The cooling time of the shocked gas depends on the 
density as 



^cool 



kT 
JMT) 



(10) 



where A( r) is the cooling fun c tion a t temperature T. Follow 
in g e.g.. iMcKee fc Ostrikel^ (|l977l ). iBlondin et all (|l998h . 



or iTruelove fc McKee (199t), A(T) can be approximated 



with piecewise power laws, with a dependence of in the 
temperature range of relevance for the shocked gas, around 
10® K. By relating the temperature to the shock velocity v, 
kT ~ rripV^ (where rup is the proton mass), and equating 
icooi to the age t as expressed in Eq. [T] one obtains that the 
trans ition radius, rmax, scales as (e.g. iBandiera fc PetrukI 
l20ld ) 



(3-2e)/(ll-6e) - (5-2e)/ (ll-6e) 



(11) 



A remnant expanding beyond this radius, at the given am- 



bient density, will enter the radiative phase and quickly fade 
from view. For a fairly large range of plausible cooling func- 
tion dependences, e.g., indices e of —1/2 to —3/2, rmax de- 
pends on density as p~^^'' to p~'^^^ . Conversely, the maxi- 
mum ambient density that will permit a Sedov-phase SNR 
of radius r is 



Pmax oc r 



(12) 



where 5 is likely in the range —7/3 to —5/2. 

The expected size distribution of remnants will now be 



dN 
dr 



dN dt 
dt dr 



,3/2 



,dP_ 
dp 



P^'^dp, (13) 



where dP/dp is the distribution of gas densities in a given 
galaxy, and pmin is the minimum density in that galaxy, on 
the size scales relevant for SNRs. Let us parametrize the 
density distribution as a power law of index /3, 

dP_ 

dp 



(14) 



with /3 > —3/2. Then, substituting in Eq. 1131 integrating, 
and expressing pmax in terms of r using Eq. 1121 gives a rem- 
nant size distribution 



dN 
dr 



^«(/3+3/2)-|-3/2 



(15) 



A uniform size distribution, dN/dr ~ const., is obtained if 
P = -6/7 (for S = -7/3) or ^ = -9/10 (for 5 = -5/2). 
In other words, the observed, roughly uniform, SNR size 
distribution implies a powerlaw distribution of densities with 
index j3 ~ —1. 

We note that, because of the dependence of rmax on 
^cooi, rmax cau depend on the metallicity of the gas, Z. As- 
suming a cooling function proportional to Z, Eq. 1111 in ad- 
dition to the dependence on p and Eo, will depend on Z 
as 



-2/{ll-6£) 



(16) 



However, in dwarf galaxies like the MCs, where Z does not 
change much from one location to another, this will be a 
small effect. For example, if e = — 1, changing Z by a factor 
5 would only change rmax by ~ 20%. 
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The existence of a cutoff in the distribution beyond 
some size is naturally expected in the above picture, given 
the existence of some minimum density pmin in the regions 
of the galaxies where SNe explode, and in view of the de- 
creeise of pmax with radius. At some radius, pmax will equal 
Pmin and the integral in Eg. 1131 will therefore become zero. 
In other words, there will be nowhere in the MCs a region 
with a density low enough to permit a SNR of that size that 
is still in its bright Sedov phase. 

A deficit of SNRs at small radii, as observed in M33 and 
the MCs, is also expected in this scenario. Before the onset 
of the Sedov stage, faster shock velocities will lead to fewer 
objects observed in the bins with the smallest radii. This 
onset happens at ages (sizes) that depend on the details of 
the ejecta structure, as w ell as the ambient density (see § 7 in 
iTruelove fc McKedl 19991 ). but for most SNRs it should occur 
around a few hundred years (a few pc), which is consistent 
with the deficits that we have observed in M33 and the MCs. 
This regime affects only a small number of objects in the 
MCs, and does not impact any of the arguments made above, 
so we will ignore it for the remainder of the paper. 

We have also ignored the deviations from the standard 
evolutionary picture that can be introduced by the shape of 
the circumstellar me dium excavated by the SN progenitors. 
I Badenes et al.l (|2007l ') showed that most Type la SNRs with 
known ages have sizes that are consistent with an interac- 
tion with a uniform ambient medium, but no such study has 
been done for CC SNRs. The fast stellar outflows expected 
from the more massive CC SN progenitors will modify the 
sizes of a fe w individual SNRs at c ertain stages of their evo- 
lution fe.g.. lDwarkadasll2005l . l2007l '). but most of the objects 
that we consider here are too large to be expanding in even 
the most extreme wind-blown cavities. As long as the bulk 
of the SNRs in the sample spend most of their lifetimes in 
the Sedov stage, this should not affect our scenario. Simi- 
larly, tlw_Ja£t_thM_somBSN^^ evolve inside superbubbles 
(e.g. iMac Low fc McCravl Il988h is naturally incorporated 
into our picture - superbubbles merely become one more of 
the factors driving the density distribution in the interstellar 
medium. Incidentally, the size distribution of superbubbles 
also relat es to the prope r ties o f the interstellar medium, as 
shown bv lOev fc Clarkj (|l997l ) for several nearby galaxies, 
including the SMC. 



5 THREE ESTIMATES OF THE GAS 
DENSITY DISTRIBUTION IN THE 
MAGELLANIC CLOUDS 

As we have seen, a uniform SNR size distribution can be 
understood as the result of Sedov expansion, combined with 
a transition to the radiative phase at an age that depends 
on the local density, provided that the density of the gas 
in the interstellar medium follows a distribution close to 
a power law with an index of —1, dP/dp ~ p^^. In this 
Section, we test this hypothesis by examining three indirect 
tracers of gas density in the Magellanic Clouds: HI column 
density; star-formation rate (SFR) based on resolved stellar 
populations; and Ha emission-line surface brightness. These 
tracers are well suited for our goals because they are valid 
over a wide range of densities, and the necessary data are 



available from public surveys that cover the whole extent of 
the Clouds, as described in detail below. 



5.1 21 cm emission-line-based HI column density 

We have taken the surface brightnes s of HI 21 cm line emis- 
sion in the MCs f r om t he maps of H im et al. 1 (|2003D and 
IStanimirovic et ahl (|l999l ). which combine single-dish Parkes 
and aperture-synthesis ATCA data to probe both small and 
large scales in the LMC and the SMC, respectively. The 
21 cm emission is optically thin, so the surface brightness 
is directly proportional to the HI column density. Since 
the LMC possesses a fairly face-on (inclination i ~ 35°, 
Ivan der Marel fc Cionil 120011 ). well-ordered HI disk, the col- 
umn density shou ld, in turn, be rou ghly proportional to the 
volume density p. iKim et al] (|2007l ) report that the HI col- 
umn density distribution of individual "clouds" of neutral 
hydrogen in the LMC follows a log-normal form, rather than 
a power law. However, their figure 13 suggests that, above a 
low cutoff of 2 X 10'^'' cm~'^, the distribution does behave as a 
power law of slope ~ — 1, over at least an order of magnitude. 
To re-examine this, in Figure[S]we show the differential dis- 
tribution of HI column densit y in the individual beam-sized 
pixels of the lKim et aH l|2003h LMC m ap (as opposed t o the 
cumulative plot for "clouds" shown in lKim et al.ll2007l ). We 
see that the HI column in the LMC does follow an index — 1 
power law fairly well, between a column of 3 x 10^" cm~^ and 
6 X lO'^^ cm^^. The observed low cutoff in the column den- 
sity is unavoidable because of the integration through the 
disk and over the beam size (every line of sight is basically 
sampling the densest regions at that point). In regions with 
low density, the H atoms might be ionize d, as in the "w arm 
ionized" phase of the interstellar medium (|Ferrierell200ll ) , so 
the tracer may become less reliable there. It is quite plau- 
sible that, in the regions where SNe actually explode, the 
underlying distribution of densities also reaches a minimum 
- as we recall, a minimum density is required in order to 
reproduce the observed upper cutoff in SNR size - but this 
does not necessarily correspond with the lower threshold in 
the HI distribution. To illustrate this, we have calculated 
th e mean HI columns in eac h of the spatial "cells" d efined 
by [Harris fc Zaritskvl (|2004l ) and [Harris fc Zaritskvl (|2009l ) 
that contain SNRs (see § 15.21 below for a description of the 
cells), which we display with the horizontal rulers in Figure 
[5] We note that, in the LMC, these cells have average col- 
umn values between 5 x 10^" cm^^ (close to, but higher than 
the low HI cutoff) and 6 x 10^^ cm"^, although most SNRs 
appear clustered around 2 x 10^^ cm^^. 

As shown in Figure O the distribution of HI col- 
umn densities in the SMC is flat over the same range, 
although the rise and fall from the plateau happen at 
the same densities as the rise and fall of the powerlaw 
in the LMC. While we do not know the reason for this, 
we speculate that it may be related to an SMC geometry 
that is elongated along our line of sight, and the integra- 
tion effect that results. The actual "depth" of the SMC, 
whether just a few kpc or as much as 20 kpc, is debated 
(iHatzidimitriou fc Hawkind 1 19891: [Harris fc Zaritskvl |2004| : 



ISubramanian fc SubramaniamI 20091 ). but is likelv at least 
a few times larger than that of the nearly face-on LMC. 
Such an integration effect would explain why the SMC SNRs 
are found at HI column densities that are a factor ~4 higher 
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Figure 5. Binned histogram of neutral H column densities in the 
LMC (red) and the SMC (blue). The dotted line overlaid on the 
LMC histogram is a power law of index —1 that intersect the 
data at a column of 5 X 10^''cm~'^. The column density values 
for the cells that contain the SNRs in both galaxies have been 
represented by tick marks on the horizontal rulers. 
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Figure 6. Binned histogram of star formation rates averaged over 
the last 35 Myr in the LMC (red, solid plot) and the SMC (blue, 
dashed plot). The dotted lines are power laws of index —1 that 
intersect the data at a SFR of 2 X lO'^ Mq yr"! cell"!. The SFR 
values for the cells that contain the SNRs in both galaxies are 
represented by tick marks on the horizontal rulers. 



than the LMC SNRs (see horizontal rulers in Figure[5]). Two 
other tracers of density, examined below, do suggest a —1 
power law in the SMC, as well as in the LMC, and do not 
display this ofTset in SNR ambient densities. 

5.2 Schmidt Law plus star formation rates from 
resolved stellar populations 

As a second way to estimate the local gas densities 
and test the powerlaw density distribution hypothesis, 
we use star formation rates (SFRs) calculated from the 
st ar-formation history (SFH ) m aps of the MCs publishe d 
by iHarrls^ fc Zaritskvf (|2004l ) and [Harris fc Zaritskvl (|2009l ). 
which are based on resolved stellar population^ The maps 
were elaborated using four-band ( U, B, V, and /) pho- 
tometry from the M agellanic Clouds Photometric Survey 
IZaritskv et al.ir2004l ). which has a limiting magnitude be- 
tween 20 and 21 in V, depending on the local degree of 
crowding in the images. In each Cloud, the data were di- 
vided into regions or "cells" with enough stars to pro- 
duce color-magnitude diagrams of suffici ent quality, which 
were then fed into the StarFISH code (|Harris fc Zaritskvl 
20011') to derive th e loca l SFH for each cell. For the LMC, 
Harris fc Zaritskvl (|2009l ) divided more than 20 million stars 
into spatial cells encompassing the central 8° x 8° of the 
galaxy (see their figure 4). In total, they produced 1376 cells 
for the LMC, most of them 12' x 12' in size, and about 50 
cells in regions of lower stellar density wit h a la rger size 
(24' X 24'). For the SMC. iHarris fc Zaritskvl (|2004l ) divided 
over 6 million stars into 351 12' x 12' cells, leaving out two 
areas that are contaminated by Galactic globular clusters in 

the foreground (see their figu re 3)^ 

iHarris fc Zaritskvl (|2004l ) and lHarris fc Zaritskvl (|2009l ') 
provide the SFH for each of these cells, which we use in Pa- 
per II to derive the delay time distribution of SNe in the 
Clouds, but for our present purposes we are only interested 

^ The complete maps are available at 

|http: //ngala. as ■ arlzona. edu/dennls/mcsurvey .html | 



in the current SFR, which we have obtained by averaging 
the star formation in each cell over the last 35 Myr. The 
ISchmidd ljl959l ) law relates star formation to gas mass col- 
umn, E, roughly as SFR oc E^''^. If gas density is distributed 
as p*^, and E oc p, then SFR will be distributed as a power 
law of index 2(/3 — l/2)/3. Coincidentally, if /3 = —1, the 
SFR distr ibution wil l also have an index —1 powerlaw dis- 
tribution. iKennicuttI (|l998l ) gives an update of the Schmidt 
law, relating star formation rate surface density, Espr, to 
gas mass column Egas, as 

/ y, \ 1.4±0.15 

Espr = (2.5±0.7)xl0-* , Mq yr-^kpc-^ 

\Me pc-V 

1—^-1 ^^^^ 
As emphasized by IKennicuttI (| 19891 ) , the Schmidt law has 

a threshold at a given mass column, below which the star- 
formation r ate falls steeply. T he exact value of this threshold 
is disputed: IKennicuttI (| 19891 ) finds Egas = 3 - 4 M© pc"^, 
which corresponds to ~ 1.4 x 10 Mq yr kpc , but 
higher and lower threshold values have als o been re- 
ported in the literature (see iBigiel et al.l |2008| . for an up- 
dated discussion of this issue). From Eq. [17] with Egas = 
3.5 M0pc~^, the threshold mass column corresponds to a 
SFR of 0.0014 Mq yr~^kpc~^ or 4.3 x 10"^ Mn yr"^ per 
12' X 12' cell in the IHarris fc Zaritskvl (|2009l ) maps of the 
LMC. 

Figure [6] shows that, above a value of ~ (2 — 4) x 
10~*Mq yr~^ cell~^, the SFR is indeed distributed as a 
power law of index —1, over more than one order of magni- 
tude in SFR, for both the LMC and the SMC. To meaning- 
fully compare the SFRs in cells with different physical areas, 
we have counted the few 24' x 24' cells in the LMC as four 
12' X 12' cells with the same SFR, and we have scaled down 
the SFRs of all the 12' x 12' cells in the SMC by a factor 
of 1.44, to account for the fact that this galaxy is 20% more 
distant than the LMC. 
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Figure 7. Binned histogram of Ha surface brightness in the LMC 
(red) and the SMC (blue). Dotted Uncs arc power laws of index —1 
that intersects the data at 2 x 10^ dR pixel" ^. The flux values for 
the cells that contain the SNRs in both galaxies are represented 
by tick marks on the horizontal rulers. 



5.3 Schmidt Law plus star formation rates from 
Ha emission 

Rather than measuring SFR via the resolved stellar popula- 
tions of the MCs, we can trace it by means of Ha emission, 
which gives us a third probe of gas density. Ha is principally 
powered by photoionisation from O-type stars, whose num- 
bers are proportional to the SFR. The Schmidt law again 
connects this SFR to the gas mass column. We have ex- 
amined the continuum-subtracted Ha emissi on maps of the 
LMC and SMC from the SHASSA survey of ICaustad et al.l 
(j200ll ) 0- Figure [3 shows that both the LMC and the SMC 
have distributions of Ha surface brightness that follow power 
laws of index ~ — 1, this time over 2 orders of magnitude for 
the SMC, and almost 3 for the LMC. As before, through the 
Schmidt law, this implies a power law with index —1 for the 
gas density distribution over 1-2 orders of magnitude in den- 
sity. Surface brightness is independent of distance for such 
nearby galaxies, and can therefore be compared directly be- 
tween the LMC and SMC. 



5.4 The Interstellar Medium of the Magellanic 
Clouds 

The density tracers that we have examined here are not 
free of biases, but they should provide a fairly accurate pic- 
ture of the density distribution in the interstellar medium 
of the MCs. Specific regions within the Clouds, like the 
LMC bar or 30 Dor, will probably have density distribu- 
tions that are different from the simple powerlaws we have 
found here; our analysis only applies to the bulk, statistical 
properties of the gas on galactic scales. Several theoretical 
studies have adressed the probabilit y distribution of den - 
sities in the in terstellar medium: sec 'Padoan et al" |l997l) . 
Passot fc Vazgucz-Scma dcni (1998), Scalo ct alTi 1998[ ). and 
Wada fc NormanI (|200li ). The consensus of these studies is 



that, despite local inhomogeneities and deviations, a com- 
bination of galaxy-wide and local mechanisms conspire to 

^ The maps are available at |http : //amundsen . swarthmore ■ edu/ 1 



produce a lognormal distribution of densities in the gas for 
a wide range of conditions. At the higher densities where 
most stars are formed and most SNe explode, the tail of 
the lognormal distribution is well repro duced by a power 
law w ith an index ~ — 1 (see figure 19 in I Wada fc NormanI 
l200'iD . This is consistent with our findings for the MCs. 



6 DISCUSSION 

The derivation in § |4] ignores many subtleties in what is 
clearly a complex problem, involving different physical pro- 
cesses and temporal scales. In this context, the fact that the 
distribution of gas densities in the MCs seems to behave like 
a power law with an index ~ — 1 does not prove that our 
proposed picture for the evolution of SNRs is correct. How- 
ever, our simple (but physically motivated) scenario appears 
to explain the available observations quite well, and is con- 
sistent with what we know about the bulk properties of the 
interstellar gas in the MCs. Furthermore, there are indirect 
ways to test that the basic scenario is sound. 

In Paper II, we use the compil ation of SNR observations 
fro m Table 1 and the S FH ma ps of lHarris fc Zaritskvl (|2004l ) 
and [Harris fc Zaritskvl (|2009l ) to derive the SN rate and de- 
lay time distribution in the Magellanic Clouds. Some key 
results from that exercise, which we describe briefly here, 
serve as plausibility arguments for our picture of SNR evo- 
lution. A basic ingredient in the derivation of SN rates is 
the "visibility time" of SNRs (i.e., the time during which 
a SNR would be visible, if it were there). We have iden- 
tifled this time with the transition to the radiative phase 
in § |4l but without assigning any specific value to it. This 
cannot be done using theoretical arguments alone, because 
many of the necessary parameters are not known with suffi- 
cient accuracy, like the distribution of SN kinetic energies or 
the normalization of the cooling function (see Eq. illjl . The 
visibility time cannot be calibrated using individual objects 
either, because the only SNRs that have reliable ages from 
light echoes or historical observations are much younger than 
the old SNRs in the Sedov stage that form the bulk of our 
sample. 

In Paper II, this conundrum is solved by imposing a con- 
dition on the delay time distribution: that the vast majority 
of stars more massive than 8Mq explode as core collapse 
SNe, with very few "silent" events that collapse directly to 
form a black hole without much ejection of material and 
hence without leaving a SNR behind. If this condition holds, 
an absolute value for the SNR visibility time can be obtained 
by equating the time-integrated SN rate in the temporal bin 
of the delay time distribution that is associated with core 
collapse SNe (in the binning used in Paper II, that is all SNe 
with delays shorter than 35 Myr) to the number of massive 
stars per total stellar mass formed. This procedure yields 
SNR visibility times that depend on the tracer used to de- 
termine the local density, and vary between 13.3 kyr (for the 
HI density tracer) and 22.5 kyr (for the SFR-based density 
tracer - see table 2 in Paper II). These ages are for regions 
where the local density corresponds to the mean value of 
the tracer, 1.5 x 10^^ cm"^ and 3.3 x 10"* Mq yr'^ ceU"\ 
respectively (see Figures [S| and [S| ; for values 10 times lower 
than the mean, the visibility times could be a factor 3 to 4 
longer (see Paper II for details). This is in good agreement 
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with the maximum ages of SNR shells inferred by SNR- 
pulsar associations in the Milky Way, which are ~ 60 kyr 
llFrail et al ] |l994l : lGaensler fc Johnstonlll995l l. 

Another key result from Paper II is the total rate of 
SNe (CO and Type la) in the Magellanic Clouds. The rate 
depends again on the density tracer used to estimate the vis- 
ibility time, and varies between 4.1 ±0.9 X 10"^ SNyr-i (for 
HI) and 2.4 ±0.4 x 10"^ SNyr"^ (using SFR as the tracer). 
These rates, and the visibility times quoted in the previous 
paragraph, are derived assuming that our SNR sample is 
complete and our evolution scenario is correct. As explained 
in Paper II, the rates, when normalized per unit stellar mass, 
are typical of dwarf irregular galaxies, of which the MCs are 
prototypes. The rates are also in loose agreement with the 
historical records, which show 2 SNe in the MCs in the last 
40 yr (i.e., ro ughly 5 x 10"^ SN yr ^, see discussion in § 5.3 
of iBadenes et a l. 2008), although this number of events is 
too small to put strong constraints on the SN rate. Adopt- 
ing the assumption that most SNRs are actually in free ex- 
pansion and not in the Sedov stage, as has been invoked 
by a number of authors to explain the observed uniform size 
distribution (see §|3}, would lead to inconsistencies. Undecel- 
erated SN ejecta reach velocities in excess of 10000 km s~^, 
so we would expect the shock velocities in freely expand- 
ing SNRs to be in the range ~ 5000 to 10000 km s~^ At 
these velocities, the cutoff radius of r ~ 30 pc would trans- 
late to maximum ages of 3 to 6 kyr. The 77 SNRs in the 
MCs would then imply a SN explosion in the MCs once 
every 40 to 80 years, on average, or even more frequently 
if our SNR sample is incomplete. Finally, if the visibility 
time were fixed at 6 kyr, for instance, this would correspond 
to a CC-SN yield (i.e., the number of CC SNe per stellar 
mass formed) of 0.022 Mq^ (see Paper II). To produce so 
many CC-SNe all stars with initial mass above 3.7 Mq in 
a standard IMF would have to explode, which is in direct 
contradiction to stellar evolution theory, an d to the semi- 
empirical initial-final mass relation for WDs ICatalan et al.l 
l2008l: ISalaris et al]|2009l: IWilhams et "aLll2009lV 



column density. Ha flux, and SFR obtained from resolved 
stellar populations. We have verified that these tracers do in- 
deed indicate a density distribution that follows a power law 
with index —1, over a wide range of densities, and that this 
agrees with our theoretical understanding of the dynamics 
of the interstellar medium. Our scenario implies that SNRs 
will remain visible for different times at different locations 
in the Magellanic Clouds, depending on the local density. 
This visibility time is a key ingredient in the calculation of 
SN rates and SN delay time distributions, which we derive 
for the Magellanic Clouds in Paper II. The absolute value 
of the visibility time cannot be determined from theoretical 
arguments alone, but we have used the results from Paper II 
to verify that the range of values we obtain (13 to 23 kyr for 
regions of mean density, depending on the tracer) is consis- 
tent with the existing limits for the maximum ages of SNRs, 
and with basic stellar evolution theory. 

It would be interesting to extend our analysis of SNR 
sizes and interstellar medium densities to other galaxies in 
the Local Group. On the one hand, this would allow us to 
test the validity of our SNR evolution scenario in differ- 
ent settings. On the other hand, the techniques developed 
here and in Paper II would yield more accurate SN rates 
and more detailed delay time distributions with larger SNR 
samples, provided the SFHs in the host galaxies could be 
determined with enough spatial and temporal resolution. A 
growing number of nearby galaxies have published SNR cat- 
alogues, but samples of the quality of the one we have assem- 
bled here for the M agellanic Clouds or the one published by 
iLong et al.l (|2010D for M33 are hard to obtain. Interstellar 
extinction and distance uncertainties plague the Milky Way 
SNRs, and the radio detection limits become comparable to 
the fluxes of the fainter Magellan ic Cloud SNRs for distance s 
beyond several hundred kpc (see lChomiuk fc Wilcotdl2009l ). 
With a judicious investment of observing time, however, rea- 
sonably complete multi-wavelength catalogues of SNRs in a 
few well-suited nearby galaxies could become available in 
the near future. 



7 SUMMARY AND CONCLUSIONS 

In this paper, we have proposed a physically motivated sce- 
nario for the evolution of SNRs in the interstellar medium of 
galaxies, and we have applied it to the sample of SNRs in the 
Magellanic Clouds. We compiled multi-wavelength observa- 
tions of the 77 known SNRs in the MCs from the existing 
literature (Table 1 and Figure [TJ. We verified that this com- 
pilation is fairly complete, and that the size distribution of 
SNRs is approximately fiat, within the allowed uncertain- 
ties, up to a cutoff at r ~ 30 pc, as discussed by other au- 
thors before. We noted that these featur es are also p resent 
in the larger SNR sample assembled by iLong et al.1 (|2010l ) 
for M33. In our model, the flat size distribution can be ex- 
plained if most SNRs are in the Sedov, decelerating, stage 
of their expansion, quickly fading below detection as soon 
as they reach the radiative stage at an age (size) that de- 
pends on the local density. Under these circumstances, a flat 
distribution of SNR sizes arises naturally if the probability 
distribution of densities in the gas follows a power law with 
index —1. To test this hypothesis, we have examined the dis- 
tributions of three different density tracers in the Clouds: HI 
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